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Elizabeth  R.C.  Marques  and  James  H.  Williams,  Jr. 
Massachusetts  Institute  o£  Technology 
Department  of  Mechanical  Engineering 
Cambridge,  Massachusetts  02139 


SUMMARY 


The  homogeneous  problem  of  stress  wave  propagation  in 
unbounded  transversely  isotropic  media  is  analyzed.  By  adopting 
plane  wave  solutions,  the  conditions  for  the  existence  of  the 
solution  are  established  in  terms  of  phase  velocities  and 
directions  of  particle  displacements.  Dispersion  relations  and 
group  velocities  are  derived  from  the  phase  velocity  expressions. 
The  deviation  angles  -  angles  between  the  normals  to  the  adopted 
plane  waves  and  the  actual  directions  of  their  propagation  -  are 
numerically  determined  for  a  specific  fiberglass^epoxy  composite. 
A  graphical  method  for  the  construction  of  the  wave  surfaces  is 
introduced,  using  magnitudes  of  phase  velocities  and  deviation 
angles.  The  simpler  results  for  the  case  of  isotropic  media  are 
shown  to  be  contained  in  the  solutions  for  the  transversely 
isotropic  media. 
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INTRODUCTION 


Most  studies  on  of  stress  waves  in  solid  elastic  materials 
refer  to  isotropic  media.  The  usual  model  adopted  to  describe  the 
geometry  of  such  waves  is  based  on  the  theory  of  sound  waves  in 
homogeneous  materials  as  first  presented  by  Lord  Rayleigh  [1]. 
The  same  concepts  of  wave  front, dispersion  relations  and 
directivity  functions  which  are  used  to  describe  sound  waves  are 
also  used  in  the  description  of  stress  waves  in  isotropic  solids, 
as  can  be  seen,  for  instance,  in  the  work  of  Miller  and  Pursey 
[2]  and  Achenbach  [3].  Because  an  isotropic  medium  has  the 
simplest  dynamic  behavior  among  all  possible  cases  of  symmetry, 
it  can  be  used  as  a  "reference"  when  anisotropic  cases  are 
studied.  Concepts  (such  as  phase  and  group  velocities  ,  velocity 
surfaces  and  wave  surfaces  )  used  to  describe  the  propagation  of 
waves  in  isotropic  media  have  also  been  applied  to  anisotropic 
media  [4-8]  . 

This  work  is  focused  on  the  problem  of  propagation  in 
anisotropic  media.  The  behavior  of  stress  waves  in  one  type  of 
anisotropic  medium,  the  transversely  isotropic  material,  is 
described  for  an  infinite  medium,  establishing  the  parallel 
between  such  materials  and  the  isotropic  case.  Also  a  graphical 
technique  for  the  construction  of  the  wave  surface  is  introduced, 
giving  a  better  understanding  of  special  geometric  features  that 
frequently  occur  during  the  propagation  in  filamentary  materials. 
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BASIC  CONCEPTS 


Definition  of  Phase  Velocity 

The  definition  of  phase  velocity  is  formulated  by  following 
the  path  of  a  small  segment  of  a  wave  front  propagating  through  a 
medium  .  The  segment  is  assumed  to  be  sufficiently  small  that  it 
can  be  approximated  by  a  plane  wave  segment.  The  concept  is  valid 
for  any  elastic  medium  , isotropic  or  anisotropic. 

The  phase  velocity  is  defined  for  periodic  or  nonperiodic 
waves  as  [9] 


V  =  d  /  t  (1) 

n 

where  d  is  the  distance  travelled  by  the  plane  segment  in  the 
direction  of  the  normal  ri  to  the  wave  front  in  time  t  (see  Figs. 
1.a  and  l.b  )  .  The  subscript  n  for  the  velocity  indicates  that 
the  phase  velocity  represents  the  speed  of  propagation  of  the 
front  segment  in  the  direction  of  the  normal  jn. 

For  periodic  waves,  the  phase  velocity  can  be  represented 
in  terms  of  either  wavelength  and  period  or  frequency  and  wave 
number  as  [8] 


v=X/T=w/k  (2) 

n 

where  X  is  the  wavelength  ,  T  is  the  period  of  the  wave,  oj  is 
the  radian  frequency  and  k  is  the  wave  number  magnitude. 

It  is"  assumed  that,  in  general,  the  plane  wave  segment 
travels  along  a  direction  that  is  different  from  the  the  wave 
front  normal  n.  If  this  direction  is  called  r  (as  for  ray  of 


propagation  ),the  angle  between  r^  and  ^  is  defined  as  the 
deviation  angle  [10].  It  is  known  that  the  deviation  angle  for 
isotropic  materials  is  zero  [9]  . 

Also,  it  can  be  shown  that  stress  waves  in  infinite  elastic 
media  are  nondispersive  [4t6]  .  As  a  consequence  ,  the  phase 
velocities  and  deviation  angles  are  independent  of  frequency  for 
any  symmetry  configuration  of  properties  of  the  medium  , 
including  the  isotropic  and  transversely  isotropic  media 
considered  here  . 
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Definition  of  Group  Velocity 


The  group  velocity  is  defined  as  the  velocity  of  propagation 
of  energy  in  a  medium  [4]* 

For  a  periodic  wave  propagating  in  a  medium, the  components 
v^^  of  the  group  velocity  in  a  cartesian  coordinate  system  x, 
y  and  2  are  [9] 

v.=daj/dk  (3) 

gi 

where  the  's  (i=x,y,z)  are  the  components  of  the  wave 

number  vector  along  the  coordinate  axes  and  w  represents  the 
radian  frequency  of  oscillation  of  the  wave  .  The  frequency  is 
usually  represented  as  a  function  of  the  wave  number  components 
and  the  elastic  properties  of  the  medium  ;  this  representation 
is  known  as  the  dispersion  relation  [11  ].  If  the  dispersion 

relation  is  known,  the  derivatives  in  eqn.(3)  can  be  easily 

calculated.  These  relations  will  be  derived  later  from  the 
expressions  of  phase  velocity  for  the  media  of  interest, 
including  isotropic  and  transversely  isotropic  materials. 

In  unbounded  elastic  isotropic  media  the  elastic  properties 
are  the  same  in  all  directions.  For  this  reason  it  can  be  shown 

that  the  magnitudes  of  the  phase  and  group  velocities  are  equal 

and  the  directions  n  and  £  are  coincident  [7]  (see  Fig.I.a). 
These  results  are  to  be  expected  since  the  system  is 
nondispersive . 

For  a  transversely  isotropic  medium,  the  magnitudes  of  phase 
velocity  and  group  velocity  are  distinct,  and  it  can  be  shown 
that  the  deviation  angle  /\  , between  n  and  r,has  the  value  [IO] 
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(4) 


cos  a  V  /v 
n  g 

where  ®  deviation  angle  A  is 

indicated  in  Fig.  1b  . 

The  physical  meaning  of  a  nonzero  deviation  angle  is  that 
energy  propagates  obliquely  with  respect  to  the  wave  front 
normal.  The  fact  that  nondispersive  systems  can  have  distinct 
values  of  phase  and  group  velocities  was  first  emphasized  by 
Lighthill  [12]. 
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PLANE  WAVE  SOLUTION  OF  THE  HOMOGENEOUS  PROBLEM 


The  solution  of  the  equations  of  motion  is  derived  for  the 
specific  problem  of  the  propagation  of  stress  waves  in 
transversely  isotropic  media.  All  the  assumptions  and  results 
also  apply  to  isotropic  media  ,as  they  are  a  subclass  of 
transversely  isotropic  media  .  In  the  following  sections  all 
results  are  first  derived  for  a  transversely  isotropic  medium  and 
then  simplified  to  an  Isotropic  medium. 


Equations  of  Motion  for  Transversely  Isotropic  Medium 

The  following  assumptions  are  made  for  the  derivation  of  the 
equations  of  motion: 

(1)  The  medium  is  homogeneous  and  has  constant  density. 

(2)  The  medium  obeys  Hooke's  law. 

The  equations  of  motion  can  be  obtained  for  any  elastic 
medium  by  the  force-dynamic  requirements  of  a  volumetric  element. 
Using  a  cartesian  coordinate  system,  oxyz,  as  the  reference 
system,  these  equations  can  be  written  as  [9] 


T  +  T  +T 

XX, X  xy,y  xz,z 

T  +  T  +T 

xy, x  yy,y  yz, z 

T  +  T  +T 

xz, x  yz,y  zz,z 


P  u 


,  tt 


P  V 


,  tt 


p  w 


,tt 


(5) 


where 


(r,s=x,y,z)  are  the  normal  (r=s)  and  shear 
x*  s 

(r#s)  stresses  with  respect  to  system  oxyz; 
u,v  and  w  are  the  displacement  components  of  a  point  in 
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the 


medium  along  the  directions  x,y  and 


z,  respectively;  and 
P  is  the  density. 

The  indexes  following  commas  denote  derivatives. 

Now,  assume  a  general  orthotropic  medium  in  which  the 
principal  axes  coincide  with  the  reference  system  oxyz.  The 
generalized  Hooke's  law  for  this  medium  is  [13] 


xx~ 

^13", z 

yy“ 

'^12'^, z* 

+ 

> 

CM 

CM 

O 

^23*, z 

zz“ 

‘=23'',y* 

S3\z 

xz“ 

yz^ 

‘^55*’'. 2 

*".y> 

xy“ 

''66(“.y 

t  ^ 

(6) 


where  the  C. .  are  the  nine  independent  elastic  constants  of  the 

i  J 

stlfness  matrix. 

Here  the  interest  is  focused  on  transversely  isotropic 

media,  which  comprise  a  subclass  of  the  orthotropic  media 

described  by  eqn.(6),  where  the  number  of  Independent  constants 

C.  is  reduced  to  five  by  the  following  constraints: 
ij 


^13~^23*  ‘'44“^55  ^66”  ^^11  “‘^12^''^^ 

In  accordance  with  eqn.(7),  the  xy  plane  is  taken  to  be  the 
isotropic  plane  for  elastic  properties. 
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The  equations  of  motion,  eqn.(5)»  can  be  written  in  terms  of 
the  displacements  u,  v  and  w  , using  eqns.(6)  and  (7)  for  the 
transversely  isotropic  medium  ,as 


C^^u  +C.„v  +C.-W  +C^^(u  +v  )+C..(u  +w  )=pu  .. 

11  ,xx  12  ,yx  13  tzx  66'  ,yy  ,xy'  44  »22  fXz'  ,tt 

C,,(u  +v  )+C.-U  +C..V  +C.-W  +C..(v  +w  )=  Pv  ..  (8) 

66  ,yx  ,xx'  12  ,xy  11  ,yy  13  ,2y  44'  ,zz  ,yz^  ,tt 


C,.(u  +w  )+C..(v  +w  )+C._(u  +v  )+C__w  s!  P  w  4.i 
44'“  ,zx  ,xx'  44  »zy  ,yy  13  ,xz  yz'  33  ,zz  ,tt 


Conditions  for  Existence  of  Plane  Wave  Solution 


Now,  assume  a  plane  wave  solution  of  the  form  [IO] 


Cu.v.w)=  fP  .P  .P  )  exof 2 if xn  +yn  +zn  -v  t)l 
\  f  f  /  '■  X  y  z  — ^ —  X  y  z  n 


(9) 


where  P  ,  P  ,  P  are  the  amplitude  components  of  a  particle 
X  y  z 

displacement  vector  along  the  coordinate  axes  x,  y  and  z, 

respectively,  corresponding  to  a  plane  wave  with  unit  normal  n; 

n  ,  n  and  n  are  the  direction  cosines  of  the  unit 
X  y  z 

normal  n  along  the  coordinate  axes  x,  y  and  z,  respectively;  v^ 
is  the  phase  velocity  defined  for  the  direction  n  ,  and  X  is 
the  associated  wavelength  with  the  selected  n  and  v^. 

If  the  assumed  solution  ,  eqn.(9),  is  substituted  into  the 

equations  of  motion,  eqn.fs),  the  following  expressions  are 
obtained  [  10] : 


o  p  p  p 

+  Cc,n  +C,,n  -  Pv  ]P  + 

c  66  y  44  2  n  ■'  x 


(CHo+C^^)n  n  P  +(C._- 
'12  66'xyy13 


•  C,,)n  n  P  =0 
44  X  z  z 
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fS3*‘=44>%"^^*‘S3‘=44^V2'’y"f‘'44<"/*''y^>*S3"z^- 


(10) 


The  condition  for  the  existence  of  the  plane  wave  solution 

can  be  expressed  by  setting  the  determinant  of  the  matrix  of  the 

coefficients  of  P  ,  P  and  P  in  eqn..(l0)  equal  to  zero  [10,14]: 

X  y  z 


[C  n  ^+C£^n  ^ 
'•  11  X  66  y 

0  0 

+C..n  -pv  ] 

44  z  n  •' 

'12  66  X  y 

(C  +C  )n  n 

13  44  X  z 

(C._+C,,)n  n 

12  66  X  y 

[C^^n  ^+C.-n  ^+C,,n  pv 

'•  66  X  11  y  44  z  n  ■' 

^^13^^44^"y^z 

=0 

(C  +C . . )n  n 

13  44  X  z 

(C  +C  )n  n 
13  44  y  z 

^)+C__n  ^-pv 
y  ^  33  z  n  ■' 

(11 ) 

Eqn.(11)  is  known  as  Christoff el ' s  equation  [5]  and  is 

written  here  specifically  for  transversely  isotropic  materials, 

considering  the  symmetry  conditions  and  relations  between  the 

elastic  constants  of  the  stiffness  matrix  given  in  eqn.(7).  The 

expression  for  the  determinant  is  a  cubic  equation  in  v^ ,  and  its 

solution  gives  three  possible  values  of  v  for  each  selected  set 

n 

of  n  ,  n  and  n  .  (  Recall  that  each  unit  normal  vector  n 

X  y  z  “ 

in  the  medium  is  uniquely  associated  with  a  set  of  components  n^, 

n  and  n  .)  Moreover,  the  elastic  properties  are  symmetric 
y  z 

with  respect  to  the  z  axis, so  the  direction  cosines  for  the 
normal  ra  can  be  expressed  in  terms  of  the  angle  between  ri  and 
the  z  axis,  as  [I4] 
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n/  +  =  sln^e  (12) 

2  2^ 

n  =  cos  0 
z 

The  advantage  of  expressing  the  direction  cosines  in  terms 
of  the  angle  0  is  that  the  expressions  for  the  velocities  can  be 
written  in  terms  of  a  single  variable. 

Expressions  for  Phase  Velocities 

By  combining  eqns.(l1)  and  (12),  the  expressions  for  the 
phase  velocities  associated  with  the  possible  plane  wave 
solutions  can  be  written  as  follows: 

-For  a  transversely  isotropic  medium 

(Vn)i  =[CggSin^0+C^^cos^0)/ P 

(Vn)ii  =[(C^^  +  C^^sln20+C^^cos20-|/^  )/2P  (13) 

where 

£  »[(C^^-C^^sln2e+(C^^-C^^)cos20]2+4(C^^+C^^)2(slnecos0)2 

-For  an  isotropic  medium 


< 

II 

[Cgg/P 

(14) 
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Where  for  the  Isotropic  medium  the  following  constraint  relations 


between  the  elastic  constants  of  the  stiffness  matrix  have  been 
applied: 

^12“^13“^23’  =  ^  1  "^1 2 

The  isotropic  medium  results  may  be  used  as  a  means  of 
checking  the  more  general  results  for  the  transversely  isotropic 

medium.  The  results  in  eqn.(14)  are  derived  from  the 

corresponding  velocity  expressions  for  the  transversely  isotropic 
medium.  Observe  that  the  values  obtained  for  the  phase  velocities 
in  eqn.(l4)  are  the  widely  known  values  of  velocities  of  shear 
waves,  (v^)j  and  longitudinal  waves,  which 

in  terms  of  elastic  properties  can  also  be  represented  by 

(G/P)^^^  and  (E/  P  respectively.  The  values  G  and  E  are 

the  shear  and  extensional  elastic  moduli , repectively. 


Particle  Displacement  Vector  Components 

By  using  the  computed  values  of  the  phase  velocities  , the 
amplitudes  of  the  particle  displacement  components  associated 
with  each  of  the  phase  velocity  vectors  can  be  found  .  This  can 
be  achieved  by  solving  eqn.(lO)  for  the  displacement  amplitude 
components  P^,  Py,  and  P^  for  each  of  the  phase  velocities. 
Because  at  each  point  there  are  three  velocities  for  each 
direction  defined  by  ii,  there  are  also  three  displacement 
amplitude  vectors  for  each  direction  n,  namely,  (P)^=(P  ,P  ,P  )  . 

corresponding  to 


(v  )-r*  (v  and  (v  ),tt»  respectively.  The  three  displacement 

vectors  are  normal  to  each  other  but  usually  neither  normal  to 
nor  coincident  with  the  direction  n  for  a  transversely  isotropic 
medium. 

In  the  homogeneous  problem,  there  is  no  unique  solution  for 
the  displacement  amplitude  components.  Each  phase  velocity 
(v^)jj  and  (^n^III  ®  corresponding  displacement  vector 

(P)j,  (P)jj  and  respectively  ,  which  represents  a 

mode  of  vibration.  Any  multiple  of  the  three  displacement 
vectors  (£)j»  ^~^III  also  a  solution. 

As  there  are  no  initial  conditions  in  the  homogeneous 
problem,  the  components  and  P^  obtained  by  solving  eqn.(lO) 

for  each  (^n^i  expressed  in  terms  of  an  arbitrary 

constant.  Thus,  the  useful  information  that  can  be  extracted 
from  the  solution  of  eqn.(lO)  for  the  displacements  consists  of 
the  direction  of  the  particle  displacement  vector  with  respect  to 
the  coordinate  axes.  The  magnitude  of  the  displacement  vectors 
(P).  is  expressed  in  terms  of  the  components  as 


(P) 


(P..  +P. 


I^p  2,1/2 

z  ' 


i=I,II,III  (16) 


Then,  the  direction  cosines  of  the  displacement  vectors  with 
respect  to  the  xyz  coordinate  system  are  [10] 


B  [  pv  ^-(n  ^+n  ^ )  Ccc-n  C  . 

n  ^  X  y  66  z  44 


12 


p 

_3L 

P 


(17) 


■)C^£.-n  2f,  1 

66  z 


B  r  Pv/-(n/.n/)C^^-n/c33.n/(C^^.C^^)^ 

'S2*<'66) 


where 


Observe  from  eqn.(l7)  that  the  sum  of  the  squares  of  the 

ratios  P  /P ,  P  /P  and  P  /P  is  equal  to  one,  meaning  that  only 
X  y  z 

two  of  the  three  components  can  be  independently  determined. 

The  isotropic  medium  is  a  particular  case  of  the 
transversely  isotropic  medium  and  where  the  particle  displacement 
directions  are  either  normal  to  or  parallel  with  the  direction 
n.  In  this  case,  the  displacement  vectors  are  aligned  with  or 
perpendicular  to  the  direction  n. 
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VALUES  OF  ELASTIC  CONSTANTS  OF  THE  STIFFNESS  MATRIX 

In  order  to  demonstrate  the  application  of  the  plane  wave 
solution  for  the  homogeneous  problem  of  wave  propagation  in 
infinite  media,  numerical  values  are  fixed  for  the  elastic 
constants  of  the  stiffness  matrix,  representing  the  two  media  of 
interest:  transversely  isotropic  and  isotropic  materials. 

Transversely  Isotropic  Material 

For  the  transversely  Isotropic  material,  there  are  five 

independent  elastic  constants  in  the  stiffnes  matrix  that  will  be 

taken  as  C..,  C^.,  C  ,  C__  and  C^^.  Note  that  C..  =  C..-2C,^  . 

11  12  13  33  44  12  11  66 

In  this  case,  the  material  is  a  unidirectional  fiberglass  epoxy 
composite  having  a  resin  content  of  36  percent  by  weight  and  a 
density  of  1850  kg/m^.  The  elastic  constants  of  the  stiffness 
matrix  were  determined  experimentally  by  ultrasonic  techniques 
from  samples  of  the  material  as  [15] 

C^^  »  10.581  X  10^  N/m^ 

C^^  =  4.679  X  10^  N/m^ 

C^^  =  40.741  X  10^  N/m^ 

C,,  =  4.422  X  10^  N/m^ 

44 

C^^  =  3.243  X  10^  N/m^ 

DO 

The  constant  C^^  was  calculated  from  the  values  of  C^^  and 
C^^  above  as  C.-=  4-098  x  10^  N/m^ .  See  Fig. 2  for  the  relation 

DO  1 2 

between  the  coordinate  axes  and  the  material's  orientation. 
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Isotropic  Material 


For  an  isotropic  material » there  are  two  Independent  elastic 
constants  in  the  stiffness  matrix  that  will  be  taken  as  and 
^12  *  representative  isotropic  material,  E  glass  was 
chosen. 

The  values  of  the  corresponding  elastic  constants  are 

C  =  82.658  X  10^  N/m^  and  C,-=  23.31  x  10^  N/m^  with  a  density 
11 

of  2540  kg/m^  [16]  .  This  material  in  fiber  form  is  one  of  the 
constituents  of  the  transversely  isotropic  material  described 
above  . 
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POLAR  DIAGRAMS  OF  PHASE  VELOCITIES 


The  values  of  the  phase  velocl'tles  obtained  from  eqns.(l3) 
and  (14) »  using  the  numerical  values  of  the  material  properties, 
can  be  presented  in  polar  diagrams.  Such  polar  diagrams  of  phase 
velocity  are  called  velocity  surfaces  [5]. 

The  velocity  surfaces  are  shown  in  Figs.  3  and  4  for  the 
isotropic  material  and  the  transversely  isotropic 
material, respectively.  One  quarter  of  each  surfaces'  intersection 
with  the  X2  plane  is  presented  .  The  surfaces  are  designated  as 
V(SH),  V(SV)  and  V(P)  corresponding  to  the  velocities  indexed  as 
I,  II  and  III,  respectively. 

Among  the  directions  defined  by  the  normal  n  used  to 
calculate  the  points  of  the  velocity  surfaces,  several  directions 
were  selected  and  presented  in.  Figs.  3  and  4.  These  directions 
are  represented  as  the  radii  from  the  origin  of  the  coordinate 
system  to  the  surfaces.  The  particle  displacement  vectors 
corresponding  to  these  directions  are  represented  by  the  arrows 
and  dots  ,  giving  an  Indication  of  the  relative  position  between 
the  particle  displacement  direction  and  the  direction  of  the 
normal  to  the  wave  front  n.  The  relative  position  of  the  particle 
displacement  direction  and  the  corresponding  direction  of  the 
normal  £  (radii  from  the  origin  )  can  be  used  to  identify  the 
modes  of  propagation.  If  these  two  directions  are  coincident,  the 
mode  is  purely  longitudinal.  If  the  two  directions  are 
perpendicular  ,the  mode  is  purely  transverse. 

Observation  of  the  V(P)  surface  in  Fig.  3  shows  that 
directions  £  and  particle  displacement  directions  are  coincident 
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for  all  directions  in  the  medium,  indicating  the  existence  of  a 
pure  longitudinal  mode  of  propagation  .  The  surfaces  V(SH)  and 
V(SV)  in  Fig, 3  show  that  directions  n  and  particle  displacement 
directions  are  normal  to  each  other  for  all  directions  in  the 
medium,  indicating  the  existence  of  two  pure  transverse  modes  of 
propagation.  So,  for  isotropic  materials  there  are  three  possible 
pure  modes  of  propagation,  one  longitudinal  and  two  transverse. 

For  the  transversely  Isotropic  material,  there  is  one 
possible  pure  transverse  mode  propagating  in  all  directions  of 
the  medium,  corresponding  to  the  V(SH)  surface  as  can  be  seen  in 
Fig. 4.  Pure  longitudinal  modes  occur  only  at  the  intersections 
of  the  V(P)  surface  with  the  coordinate  axes  (x,  y  and  z)  where 
the  particle  displacement  direction  coincides  with  the  direction 
n.  For  the  V(SV)  surface  there  are  three  directions  where  pure 
transverse  modes  occur,  namely,  the  directions  of  the  coordinate 
axes  X,  y  and  z.  (In  Fig. 4  only  the  Intersections  with  the  x  and 
z  axes  are  shown) . 

According  to  the  convention  adopted  for  the  designation  of 
the  surfaces.  Table  1  shows  the  modes  of  vibration  for  points 
along  the  directions  of  the  coordinate  axes  and  the  surfaces 
that  contain  the  indicated  modes  .  Other  designations  for  the 
surfaces  can  be  found  in  the  literature  [10,14»17],  as,  for 
instance,  "transverse"  for  the  SH  waves,  "quasi-transverse"  for 
the  SV  waves,  and  "quasi-longitudinal"  for  the  P  waves  . 
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DISPERSION  RELATIONS  AND  GROUP  VELOCITIES 


After  the  phase  velocities  for  each  direction  n  are 
obtained, the  group  velocities  can  be  determined  by  the  use  of 
the  dispersion  relation.  Combining  the  definition  of  phase 
velocity,  eqn.(2),  and  the  expressions  obtained  for  the 
velocities,  eqns.(l3)  and  (14)»  the  dispersion  relations  can  be 
written  as  follows: 

-For  a  transversely  isotropic  medium 

(  CO  )j  . 

(  (Ji  =  {[C  .  .k^  +  C  .  .  (k  ^  +  k  ^)+C__k  e  k4]/2  P  (18) 

II  ^  44  11 '  X  y  33  z 

(  oj  )___  =  {[C  .,k^+C.  .  (k  ^+k  ^)+C__k  ^+\/e  k^J/a 

III  ‘ 44  11'  X  y  33  2  '  ■"  ’ 

-For  an  isotropic  medium 

(  CO  )i  .  [Cg5(k/*ky^k/)/0]^''2 

(  CO  )ii  .  {Cgg(k/*k^^k/)/P  (19) 

(  CO  =  tC,,(k/.k/*k/)/P]’/2 

where  the  square  of  the  wave  number  magnitude  is  represented  by 

k^=  (k  ^+k  ^+k  (20) 

X  y  z  '  '  ' 

and  the  relations  between  the  components  of  Ic,  namely, 

k  k  ^  =  k^sin^e,  k  ^  =  k^cos^e  (21) 

X  y  z 

were  used  (as  in  eqn.(l2)  for  n) . 
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From  the  dispersion  relations  in  eqns.(l8)  and  (19)  and  the 
definition  in  eqn.(3)»  the  values  of  the  group  velocity 
components  along  the  coordinate  axes  can  be  determined  .  The 
expressions  for  the  group  velocity  components  in  the  xz  plane  are 
as  follows  : 

-For  a  transversely  isotropic  medium 

(v  ),.  =  n  [(C..n  ^+C.,n  P 

'  gx*^!  X  66  66  X  44  z  ' 

(=44-S3>"/’"^"x"z"(=13*=44>"> 

(22) 

<S4-S3>''x^^^^'’z''x^(S3"=44>^' 

'=44-S3>''/^"2''x”x"(=13"'=44>"' 

"  ■''^\<'=44-S3)f'Sl-‘=44'"x^-‘ 
(=44-S3^''x"J"2"z''/'=13""44>"' 

»here  B.{  ^ ]/2  p 
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^  =  ^  +  "  +  ^  +  c__n  ]/2P 

44  X  z  '  11x  33  z' 


-For  an  isotropic  medium 


‘V^II  '  “*'66  t'^66  " 

('gz^II  =  ”z'=66  ['=66  ^  (»/+>'2^)r''^^ 

(v  )ttt=  ^  +  ^ 

'  gX'^III  X  11  '■  11  '  X  Z  '  ■' 

(v  )ttt=  t*-n  P  (“  ^  +  “ 

'  gz  III  Z  11  '■  11  '  X  z  '  ■' 


(23) 


Observe  that  the  magnitudes  of  the  group  velocities  can  be 
obtained  from  the  corresponding  components  in  the  xz  plane  as 


(vj 


rr'  T  TT 


=  rr 


0  ^1' 


(24) 


For  an  isotropic  medium,  the 
velocities  are  identical  to  those  of 
defined  in  eqns.(l4)  • 


magnitudes 
the  phase 


of  the  group 
velocities  as 
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CALCULATION  OF  DEVIATION  ANGLES 


The  purpose  of  calculating  the  deviation  angles  is  to 
determine  the  actual  direction  r  of  propagation  of  plane  wave 
segments  having  normal  ti.  The  deviation  angle  can  be  expressed 
in  terms  of  the  phase  and  group  velocity  magnitudes  as  given  in 

eqn.(4) • 

Because  z  is  the  symmetry  axis  in  the  transversely 
isotropic  medium  under  consideration  ,  the  directions  of  phase 
velocity  and  group  velocity  (ri  and  r,  respectively)  are  expressed 
in  terms  of  their  angles  with  respect  to  this  axis  . 

For  the  phase  velocity  direction  ,  the  angle  0  with 
respect  to  the  z  axis  (see  Fig. 1b  )  is  defined  by  the  direction 
cosines  of  n  as  given  in  eqn.(l2)  as 


( tan  9 

r/  2  2x1/2,  T 

ttt=  L(n  +n  )  /  n  Jt 

,11, III 

(25) 

The 

angle 

0'  of 

the  direction  r  (for 

which  the 

group 

velocity 

is  defined) 

with  respect  to  the 

z  axis 

can  be 

expresse4  in  terms  of  the  group  velocity  components  as 


gy 


^gx^I,II,III 


The  indexes  I, II  and  III  in  eqns.(25)  and  (26)  denote  that 
the  same  equations  apply  to  each  of  the  velocity  sets.  Moreover, 
due  to  symmetry,  n  and  v  may  be  taken  as  zero  ,and  the 

y  g  y 

calculation  can  be  done  for  the  xz  plane  . 

The  general  expression  for  the  deviation  angle  is  then 


A  =  9 


(27) 
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The  specific  algebraic  expressions,  expressed  in  terms  of 
elastic  constants  for  the  deviation  angles  of  transversely 
isotropic  media,  are  not  reproduced  here.  For  isotropic  media 
the  deviation  angles  are  zero.  It  is  interesting  to  note  that 
algebraic  manipulations  show  that  the  deviation  angle  is 
independent  of  the  frequency.  And,  since  the  phase  velocity  is 
also  independent  of  frequency  (see  eqn3.(13)  and  (14)) »  so  is 
the  group  velocity,  which  is  another  way  of  expressing  the 
nondispersive  character  of  the  elastic  media  under  study. 

The  calculated  deviation  angles  for  the  transversely 
isotropic  medium  under  consideration  are  shown  in  Figs.  5a  ,5b 
and  5c  ,  for  the  SH,SV  and  P  modes,  respectively. 


22 


GRAPHIC  CONSTRUCTION  OF  WAVE  SURFACES 


In  this  section  the  graphic  construction  of  the  wave  front 
will  be  described,  focusing  on  the  transversely  isotropic 
material.  The  quantities  needed  for  this  construction  are  the 
magnitudes  of  the  phase  velocities  ,the  directions  ^  for  which 
these  velocities  are  defined  and  the  deviation  angles. 

In  other  words  ,  the  locus  of  the  points  of  equal  phase  in 
geometric  space  for  unit  time  will  be  determined,  supposing  that 
plane  waves  passed  through  the  origin  at  time  t^O  in  all  possible 
directions  defined  by  their  normals  n.  The  construction  is 
described  for  the  xz  plane  but  it  is  valid  for  any  plane 
containing  the  symmetry  axis  z. 

The  following  steps  must  be  followed  for  the  construction  (see 
Fig. (6)); 

(1)  Choose  a  direction  n  (n  ,  n  ) 

X  z 

(2)  Draw  a  line  from  the  origin  of  the  coordinate  system 

in  the  direction  ri  and  having  a  scaled  length  equal  to 

V  .  Since  time  is  taken  as  unity,  the  segment  represents  a 
n 

distance . 

(3)  From  the  line  with  direction  n,  mark  the  deviation 
angle  ,  if  positive  below  ,if  negative  above  the  line 
having  the  direction  n.  The  direction  so  determined  is  r 
which  is  the  direction  of  the  the  group  velocity  v  . 

s 

(4)  From  the  tip  of  the  line  of  length  v^ ,  draw  a  line 
perpendicular  to  it. 

(5)  The  intersection  of  the  perpendicular  line  with  the 
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direction  _r  determines  a  point  P  on  the  wave  surface. 
Thus,  point  P  represents  the  location  at  time  t=1  sec  of  a 
point  of  the  wave  front  which  at  time  t=0  passed  through 
the  origin  with  velocity  v^  in  the  direction  of  n. 

If  the  construction  described  is  repeated  for  each  set  of 
calculated  phase  velocities  assuming  all  possible  directions  n 
of  the  wave  normal  in  the  medium  ,  the  wave  surfaces  are 
obtained.  The  construction  steps  of  the  wave  surfaces  W(SV)  and 
W(P)  corresponding  to  velocity  surfaces  V(SV)  and  V(P) 
respectively,  are  shown  in  Figs.  7  and  8  for  some  directions  in 
the  representative  transversely  Isotropic  medium.  In  these 
figures,  the  primed  points  are  points  of  the  extremities  of 
the  segments  of  length  v^  as  described  in  step  2).  The  umprimed 
points  are  the  corresponding  points  on  the  wave  surface. 

The  construction  of  the  wave  surface  W(SH)  corresponding  to 
velocities  V(SH)  is  not  shown  since  the  deviation  angles  are 
small  and  the  final  shape  of  W(SH)  is  very  similar  to  the  shape 
of  V(SH).  The  three  wave  surfaces  obtained  for  the 

representative  transversely  isotropic  medium  are  shown  together 
(positive  x-z  quadrant  only)  in  Fig.  9  • 

For  the  isotropic  medium  no  construction  is  made  since  the 
deviation  angles  are  zero.  For  such  a  case,  the  wave  surfaces  are 
identical  to  the  corresponding  velocity  surfaces. 

It  is  important  to  emphasize  that  the  wave  surfaces  are  the 
actual  geometric  positions  of  wave  fronts  for  the  three  possible 
modes  of  propagation.  Because  plane  wave  solutions  are  assumed  to 
generate  these  surfaces  , the  wave  surfaces  represent  the  envelope 
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of  all  possible  plane  waves  emanating  from  the  origin 
coordinate  system  in  all  possible  directions,  as  if 
source  existed  at  the  origin. 


of  the 
such  a 
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SPECIAL  FEATURES  OBSERVED  IN  THE  WAVE  SURFACE  W(SV 


As  it  can  be  observed  in  Fig.  7,  for  the  representative 
transversely  isotropic  medium  there  is  a  folding  of  the  W(SV) 
wave  surface  due  to  an  inversion  of  the  sign  of  the  deviation 
angle 

To  explain  the  geometric  shape  of  the  wave  surface  W(SV), 
assume  that  the  wave  front  is  constituted  of  many  small  segments 
of  plane  waves  having  the  same  length  and  passing  through  the 
origin  at  time  t=0.  Assume  that  for  an  infinitesimal  interval  of 
time  these  segments  of  plane  wave  are  still  very  close  to  the 
origin  and  form  a  polygonal  line  that  can  be  Inscribed  in  an  arc 
of  circumference.  Then  immediately  thereafter,  the  segments  start 
to  deviate  from  their  normals,  following  the  directions  £  defined 
by  the  deviation  angles  .If  it  is  assumed  that  the  segments' 
extremities  remain  connected  physically  ,  the  wave  front  must  be 
continuous  in  time  as  these  segments  will  be  stretched  as  they 
travel  away  from  the  origin. 

The  deviation  angles.  Fig.  5b,  can  be  associated  with  a 
degree  of  stretching  of  the  original  straight  segments  .  The 
larger  the  deviation  angle  ,the  larger  the  degree  of  stretching 
undergone  by  the  segment.  Recall  that  the  angle  9  represents 
the  angle  of  the  normal  n  (to  each  segment)  with  respect  to  the  z 
axis . 

Observe  that  the  curve  of  has  two  maxima,  one 

positive  at  6  =  33°  and  one  negative  at  0  =  77.5°>  and  passes 
through  zero  at  9  =  60°.  If  the  curve  is  divided  in  four 
regions  as  shown  In  Fig.  10b  ,it  can  be  observed  that  in  region  1 
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the  segments  are  stretched  increasingly  until  the  first  maximum 
at  e  =  33°.  Then  for  region  2  the  segments  are  less  and  less 
stretched  until  the  deviation  angle  reaches  the  value 
zero,  represented  by  point  C  in  Fig,  5»a.  The  segment 
corresponding  to  the  zero  deviation  is  not  stretched.  For  region 
3,  again  the  segments  are  stretched  from  zero  to  the  second 
maximum  at  9  =  77.5°  and  finally  in  region  4  they  are  less  and 
and  less  stretched  until  reaches  zero  for  the  segment  whose 
normal  is  located  at  90°  with  respect  to  the  z  axis. 

In  the  model  of  the  polygonal  line  constituted  of  straight 
segments,  if  points  A  and  E,  Fig,  10, a,  are  the  intersections  of 
this  line  with  the  coordinate  axes,  for  an  infinitesimal  time 
after  zero  time  , these  points  remain  on  the  coordinate  axes  since 
the  corresponding  deviation  angles  are  zero.  Later  ,at  a  time 
t  equals  unity,  these  points  will  be  at  positions  A'  and  E', 
respectively.  Considering  that  at  the  same  time  all  segments  are 
supposed  to  remain  connected  ,the  segments  originally  in  regions 
1  and  2  (Fig.  10. b)  will  undergo  an  effect  of  pulling  away  from 
the  z  axis;  the  segments  originally  in  regions  3  and  4  on  the 
other  hand  ,  will  undergo  an  effect  of  pulling  away  from  the  x 
axis,  moving  closer  toward  the  z  axis. 

It  has  been  determined  that  the  segments  which  are  the  tips 
of  the  cuspidal  edges  B',D'  of  the  wave  surface  W(SV)  have 
normals  n  at  angles  of  54.02°  and  75.06°,  respectively,  with 
respect  to  the  z  axis.  For  time  t  equals  unity,  the  location  of 
these  two  segments  is  defined  by  the  lines  passing  through  the 
origin  and  the  points  B'  and  D',  respectively,  at  61,86°  and 
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41.93°  with  respect  to  the  z  axis  (see  Fig.  10a  ). 

The  originally  equal  segments  that  at  time  t  equals  unity 
are  located  inside  the  region  limited  by  the  lines  OB ' and  OD '  can 
be  associated  with  the  amount  of  energy  that  is  concentrated  in 
this  region.  If  it  is  assumed  that  each  segment  that  passes 
through  the  origin  at  time  t=0  carries  the  same  amount  of  energy, 
the  percentage  of  the  total  number  of  the  segments  contained  in 
the  region  of  interest  at  time  t  equals  unity  is  also  the 
percentage  of  energy  contained  in  the  region  at  any  other  time 
t.  Using  this  concept,  it  was  found  numerically  that  64.4^  of 
the  total  energy  propagates  through  the  region  located  between 
the  lines  OB'  and  OD '  where  the  included  angle  B'OD*  is  20.93°. 
This  calculation  was  made  for  a  total  number  of  180  segments. 

Concerning  the  geometry  of  propagation  of  the  entire  wave 
front,  it  can  be  said  that  if  an  observer  stands  at  a  position 
with  respect  to  the  coordinate  axes  inside  the  region  defined  by 
lines  OB'  and  OD '  (for  instance,  at  position  P  in  Fig.  10. a  ), 
such  an  observer  will  see  three  portions  of  the  same  front  W(SV) 
passing  through  him  at  different  times  . 

It  must  be  stressed  that  the  special  feature  described  above 
is  valid  for  the  fiberglass  epoxy  composite  used  in  this 
numerical  example.  The  existence  of  this  special  behavior  is  not 
a  common  characteristic  of  all  transversely  isotropic  media  and 
depends  on  the  numerical  values  of  the  elastic  constants  ,  so 
such  features  must  be  analyzed  and  determined  for  each  specific 
material . 
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CONCLUSIONS 


The  problem  of  the  geometry  of  propagation  of  waves  in 
transversely  isotropic  media  was  studied.  The  solution  of  the 
homogeneous  problem  of  motion  for  a  fiberglass  epoxy  composite 
was  analyzed  and  the  numerical  values  of  the  results  were 
presented. 

The  solution  of  the  homogeneous  problem  of  motion  in  an 
infinite  medium  provided  a  description  of  the  behavior  of  plane 
stress  waves  in  terms  of  velocities  and  paths  of  propagation 
Using  these  concepts,  a  graphical  method  for  the  construction  of 
the  wave  fronts  was  introduced  . 

As  described  above, the  following  conclusions  can  be  stated: 

(1)  The  graphical  technique  is  found  useful  in  the  construction 
and  the  physical  interpretation  of  the  wave  surfaces  .  The 
technique  was  applied  to  transversely  isotropic  materials  but 
the  concept  can  be  extended  to  other  cases  having  anisotropic 
symmetry . 

(2)  The  wave  surfaces  corresponding  to  a  fiberglass  epoxy 
composite  material  were  constructed  .Special  features  were  found 
in  the  W(SV)  surface. 

(3)  It  was  found  that  energy  carried  by  the  SV  waves  in  the 
fiberglass  epoxy  composite  material  is  concentrated  along  certain 
preferential  directions.  This  phenomenon  can  be  explained  in 
terms  of  the  spreading  geometry  of  plane  wave  segments  . 

The  conclusions  above  suggest  some  possible  applications  and 
recommendations  for  ultrasonic  NDE  of  transversely  isotropic 
materials  .  For  instance,  the  special  features  of  the  wave 
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surface  W(SV)  for  the  fiberglass  epoxy  material  can  be  used  as 
an  auxiliary  criterion  for  the  placement  of  transducers  when 
experiments  are  designed  .  For  example,  the  locations  where  the 
energy  is  focused  (angle  B’OD'  in  Fig. 10  )  can  be  selected  as  the 
receiving  region  since  the  signals  within  this  region  are 
expected  to  be  stronger  due  to  larger  displacement  amplitudes. 
Also,  certain  locations  in  the  medium  can  be  selected  as 
receiving  points  where  up  to  five  different  wave  front  arrivals 
may  occur,  when  the  three  modes  of  propagation  are  considered. 

The  knowledge  of  the  geometry  of  the  spreading  can  also  be 
used  for  crack  detection  since  the  number  of  expected  arrivals 
can  be  computed  for  each  direction  in  the  medium.  Thus,  the 
absence  of  the  arrival  of  one  of  the  expected  portions  of  a  wave 
front  might  indicate  that  an  obstacle  was  in  the  path  of  that 
portion  of  the  front.  Such  an  obstacle  might  indicate  the 
presence  of  a  flaw. 
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Table  1.  Designation  of  Velocity  Surfaces  and  Modes  of  Vibration 
at  Intersections  with  Coordinate  Axes. 
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Position  of  material  axes  with  respect  to  coordinate 
axes  for  transversely  isotropic  medium. 
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Fig .  4 


Phase  velocity  polar  diagram  for  fiberglass  epoxy 
composites  for  positive  x— z  quadrant.  (Arrows 
represent  direction  of  particle  displacements.) 
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Fig. 7'  Construction  of  wave  surface  W(SV)  from 
corresponding  velocity  surface  V(SV)  and  deviation 
angles  (^)jj  for  fiberglass  epoxy  composite  , 
for  positive  x-z  quadrant.  (Points  1  through  13 
are  on  wave  surface  :  points  1' through  13'  are  on 
velocity  surface.) 
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Construction  of  wave  surface  W(P)  from 
corresponding  velocity  surface  V(P)  and  deviation 
angles  (AJttt  fiberglass  epoxy  composite  , 

for  positive^  x-z  quadrant.  (Points  1  through  10 
are  on  wave  surface  :  points  1' through  10'  are  on 
velocity  surface.) 


Distance  traveled  by  wavefront  in  z  direction 
in  I  second,  z(m) 


Deviation  angle 

(A)„  (°) 


a)  Wavefront  assumed  shapes  for  time  t  =  e  (ABCDE) 
and  t=unity  (A'B'C'D'E* ) . 


0  20  40  <0  *0 


Angle  of  n  with  respect  to  z  axis,  ©  (°) 

b),  Deviation  angles  (^)__  showing  regions  and 
characteristic  points?^ 


Fig. 10  Model  of  geometrical  spreading  for  SV  wave. 
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